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Fractional excitation and its boundstate for the S=l/2 antiferromagnetic zigzag 

spin chain in a magnetic field 
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^Department of Mechanical Engineering, Fukm University of Technology, 6-1 Gakuen 3, Fukui 910-8505, Japan. 

We demonstrate that domain-wall type excitations carrying the fractional magnetization 
5^ — ±1/3 and their bound state describe the low-energy physics around the 1/3 plateau 
of the antiferromagnetic zigzag spin chain in the strongly frustrated region. In particular we 
discuss the relevance of such domain-wall excitations to the characteristic properties of the 
magnetization curves: the 1/3 plateau accompanying the spontaneous symmetry breaking of 
the translation, the cusp singularities above and/or below the plateau, and the even-odd effect 
in the magnetization curve. 
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1. introduction 

In order to understand the role of the frustration in 
low-dimensional quantum spin systems, the S = 1/2 an- 
tiferromagnetic zigzag spin chain has been one of the 
most fundamental quantum spin models. The Hamil- 
tonian of the model is given by 



'H — ^[Ji^i • Si+i -f J2Si ■ Si+2] — H'^^S^, 



(1) 



where S is the S = 1/2 spin operator, H is the mag- 
netic field, and Ji and J2 denote the nearest and next 
nearest neighbor couplings respectively. We also intro- 
duce a notation a — J2/J1 for simplicity. A remarkable 
point of the Hamiltonian (1) is that, although its form 
is very simple, it captures some characteristic aspects 
of physics induced by the frustration; At zero magnetic 
field, the well-known transition between the spin fluid 
and dimer phases occurs at a ~ 0.2411.^ Moreover, the 
ground state is exactly represented as the dimer state 
at a = 0.5, which is called Majumdar- Ghosh point. ^ In 
addition, we should note that the zigzag chain is real- 
ized as SrCu02,i° Cu(ampy)Br2," (N2H5)CuCl3i2 and 
F2PIMNH.13 

Recent evolution of the zigzag chain physics is that the 
magnetic phase diagram is actually obtained in Ref . [^'^] . 
Of course, a lot of efforts had been made so far,^^"^^ but 
it was a difficult problem to illustrate reliable magnetiza- 
tion curves in the strongly frustrated region. As is shown 
in FIG.l(see also FIG. 2), the obtained phase diagram re- 
veals that a variety of exotic behaviors is cooperatively 
induced by the frustration and the magnetic field. For 
0.56 < a < 1.25,^° the magnetization plateau appears 
at 1/3 of the full moment, accompanying the sponta- 
neous breaking of the translational symmetry by period 
three. Moreover, the cusp singularities in the magnetiza- 
tion curve show quite interesting behaviors around the 
1/3 plateau; As a is increased, the high-field cusp merges 
into the 1/3 plateau at a ~ 0.82. Also the low-field cusp 
merges into the 1/3 plateau at a ~ 0.7, but it appears 
again when a > 0.7. Here we note that the behaviors 



of the cusps near the saturation field or the lower crit- 
ical field are explained well by the shape change of the 
spin wave or spinon dispersion curves, from which the 
two component Tomonaga-Luttinger(TL) liquid can be 
expected for TL2 in FIG.l."'^*' In addition, we find an 
interesting even-odd effect in the magnetization curve of 
a finite size system for a > 0.7. 




J2/Ji(a) 



Ji/J2(a ) 



Fig. 1. Magnetic phase diagram of the zigzag spin chain. D: 
dimer gapped phase, P:l/3 plateau, TLl: one component TL 
liquid, TL2: two component TL liquid, and EO: even-odd be- 
havior blanch. The open circles denote the position of the cusp 
singularities. The solid circles are the boundary of the dimer 
gapped phase. The triangles mean the upper and lower edge of 
the 1/3 plateau, and the crosses indicate the upper edge of the 
even-odd effect in the magnetization curve. The saturation field 
is normalized to be unity. 



A series of the above intrinsic structural changes of the 
magnetization curves indicates that the zigzag chain still 
contains quite rich physics. In order to discuss the mi- 
croscopic mechanism behind such exotic magnetic phase 
diagram, we would like to emphasize the following fea- 
tures of the zigzag chain: 

• Frustrating interaction: Usually, the frustration ef- 
fect is discussed in the context of the triangu- 
lar structure of the lattice. However, here we 
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want to stress the competition between the sin- 
gle chain structure (J2 S> Ji) and double chain 
structure(J2 <C Ji) of the zigzag lattice, which pro- 
vides some essential insights for the frustration effect 
in the zigzag chain. 

• Translational invariance: the Hamiltonian (1) con- 
nects the above two limits continuously without loss 
of the translational invariance. The parameter a in- 
terpolates the single and double chain natures con- 
tinuously. 

• The 1/3 plateau is extending up to the Ising 
anisotropic limit, where the ttJ, order is realized. 
This means that the tti order in the Ising limit can 
be a good starting point to discuss the isotropic case. 

On the basis of the above noted features of the zigzag 
chain, in the following, we want to review that the 
domain- wall(DW) excitations carrying fractional value of 
successfully explain the characteristic behaviors of the 
magnetization curves. In §2 and §3 we briefly overview 
the important features of the magnetization curves and 
the 1/3 plateau state for the zigzag chain of some typical 
parameters. In §4, we introduce the zigzag XXZ chain, 
which captures more clearly the above-noted three key 
points, and then explain the essence of the fractional- 5'^ 
DW excitations. In §5, we describe the implications of 
the DW picture for the magnetization curves and finally 
summarize the conclusions in §6 

2. magnetization curves 

In FIG. 2, we show the magnetization curves for a = 
0.6, 0.7, 0.8, 1.0 and = 0.8, which were obtained with 
the DMRG^e for 192 spins. In the figures, the scale of 
the magnetic field is normalized by the saturation field 
Hg and m is the magnetization per spins. We also de- 
fine M as the total of the system for later conve- 
nience. First of all, we can see that the 1/3 plateau 
appears for 0.56 < a < 1.25, accompanying the spon- 
taneous symmetry breaking of the translational invari- 
ance. Here, we note that an anomalous step in the 1/3 
plateau for a = 0.7 and 0.8 is due to the open bomid- 
ary effect. In addition, the cusp singularities can be seen 
in low and high field regions. A more interesting point 
is that those cusps behave interestingly around the 1/3 
plateau; At a = 0.7, the low field cusp merges into the 
1/3 plateau and the high field cusp also merges into the 
1/3 plateau near a ~ 0.8. Moreover, for q: > 0.7, we 
can find that anomalous even-odd oscillation with re- 
spect to the magnetization M. For a > 0.82 the even- 
odd behavior is extended above the 1/3 plateau. As a. is 
increased further, the shape of the magnetization curve 
approaches that of the (double) Heisenberg chain. In the 
figure of ~ 0.8, the magnetization curve becomes 
continuous around m = 1/3, implying that there is no 
1/3 plateau. Finally in —>■ limit, the magnetiza- 
tion curve approaches to that of the Heisenberg chain. 
Here we should make a comment on the even-odd oscil- 
lation of the magnetization curve. The overhanged steps 
corresponding to M=:odd must be skipped in the true 
magnetization curve of the finite size system. However 
we have shown them, since they are reflecting an inter- 
esting aspect of the double-chains nature of the system. 



Author Name 

The properties of the cusp singularities near the satu- 
ration fleld or the lower critical field are explained well by 
the shape change of the spin wave or spinon dispersion 
curves. ^^'^^ However, these previous picture of the cusps 
is not satisfactory for a series of the subtle behaviors of 
the cusps around the 1/3 plateau. 




Fig. 2. Magnetization curves for a = 0.6, 0.7, 0.8, 1.0 and a ^ = 
0.8. 



3. 1/3 plateau state 

Let us discuss the nature of the 1/3 plateau, which is 
the starting point of analyzing the microscopic origin of 
the magnetic phase diagram. According to the quanti- 
zation condition of the magnetization plateau, the 1/3 
plateau state satisfies g(l/2 — m) =integer with m = 1/6, 
where q is the period of the plateau state. This implies 
that the translational symmetry of the Hamiltonian (1) 
must be spontaneously broken by g = 3, 6, • • • . The spin 
profile of the plateau along the chain direction is shown in 
FIG. 3, where the "tti" typo long-range order is realized. 
Moreover, the boundary effect decays rapidly, which sup- 
ports that the 1/3 plateau state has the excitation gap. 

Since the zigzag chain is a system having both of the 
geometrical frustration and the quantum fiuctuation, one 
may wonder why such a classical ttJ, order is still stand- 
ing against the strong quantum fluctuation and the frus- 
tration. This may be a very naive but important ques- 
tion to understand the zigzag chain. An important point 
is that the 1 /3 plateau phase is extending up to the Ising 
anisotropic limit, for which the tti array was proven in 
Ref. [^*] . The triangular structure of the lattice stabilizes 
the 1 /3 plateau of the tti order in wide parameter region 
in the classical limit. 

An qualitative explanation for the stability of the tti 
order in the quantum case may be possible by consider- 
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Ji=1.0, J2=0.8, N=192, M=32, E=-57.433 spin array 5^ energy 
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rpg^^jg rpjjg energy excitation on the 1/3 plateau state, 
which is described by the domain wall made from three types 
I of the Neel ordered spin aligns. These domain-walls can be re- 

96 j 192 garded as quasi-particles carrying fractional 5^. 



Fig. 3. The distribution of the local spin moment at M=32(l/3 
plateau) for a = 0.8^^ 

ing the role of the XY term. Since the XY term tries to 
change a tJ, pair on a bond into a singlet, the influence 
of the quantum fluctuation is the most signiflcant at the 
zero magnetic field(subspace of M = 0), where the num- 
ber of the ti pair is maximum. At the 1/3 plateau, how- 
ever, the number of the ][ pair is diluted by the ffj so 
that the effect of the quantum fluctuation should be re- 
duced certainly. Here we should recall that, in the M = 
subspace of the usual XXZ chain, the classical Neel order 
survives up to the isotropic limit. Thus it may be natu- 
ral that the tti order at the 1/3 plateau of the zigzag 
XXZ chain survives against the quantum effect up to the 
XY-Hke region. 

4. domain wall excitations 

In order to discuss the nature of the excitations around 
the 1/3 plateau, we introduce the Ising-likc XXZ chain 
with the nearest neighbor (NN) coupling Ji and the next- 
nearest neighbor(NNN) one J2 in a magnetic field H: 

H = H[) + Tizccman, (2) 

where 7io = TYnn + T^nnn is the Hamiltonian of the 
system part and Wzeeman = — ^^X^i^^i '^i the Zeeman 
term. The Hamiltonian of the NN and NNN interactions 
are defined as 

6N 6JV 

Wnn = Jl''^hi^^-^-l, "Hnnn = J2 fej,»-i-2, (3) 

i=l i=l 

h,, = eiS:S^ + S!sy) + StS] (4) 

where e denotes the anisotropy of the XY term. We as- 
sume the system of 6N sites so that M = N corresponds 
to the 1/3 plateau. The excitation of the usual XXZ 
chain at the zero field is often called "spinon" , which is a 
domain- wall type excitation carrying = 1/2.^°'^^ For 
the present 1/3 plateau of the zigzag chain, we show that 
the DW excitations carrying fractional = ±1/3 and 
their bound state can successfully explain the intrinsic 
aspects of the low energy physics. 

Let us first consider the classification of the excitations 
on the 1/3 plateau in the Ising limit. At the 1/3 plateau, 
one of three types of tti orders is selected, according to 
the spontaneous symmetry breaking of the translation. 
Then the excitations on such one- dimensional order arc 
described well by the DWs constructed by combinations 
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Fig. 4. 5^ = 1/3 DW excitation in the Ising limit. 

of 3 types of the ||| arrays. In Table 1, we list up 6 
possible DWs, which carry fractional value of the mag- 
netization. 

In order to see the origin of the fractional value of the 
magnetization, let us start with the single spin flipped 
state in the 1/3 plateau, which is shown in FIG. 4. The 
total magnetization of this state is increased by 1 from 
the 1/3 plateau. Using the perturbation of the XY term, 
we can systematically construct the low-energy excita- 
tions; by exchanging the nearest neighbor fj spin pair, 
we can see that the five up spin cluster is finally decom- 
posed into three ttt clusters, which are identical to the 
S"^ = 1/3 DW in Table I. Since these three DWs are 
equivalent, we can see that each DW carries = 1/3. 
Moreover, the tTTT cluster in FIG. 4 is decomposed into 
two TTT DWs, implying that the TTTT DW in Table I is 
regarded as a bound state of two = 1/3 DWs. 

We analyze the one-body problem of the = 1/3 
DW in the infinite length chain within the first order of 
e. Although the higher order terms are required for a 
quantitative analysis of the problem, we can capture the 
intrinsic property of the DW excitation within the first 
order theory. Let us label the position of the = 1/3 
DW by X which is defined as a site of the center of three 
up spins. Then the matrix element of Wo in the low- 
energy region is obtained as 

Ho\x) = ^-^^\x) + £^ (la; - 3) + |x + 3)) . (5) 

Here an important point is that the NNN term generates 
only higher energy configurations accompanying the en- 
ergy rise of order Ji or J2. Thus we can obtain the low- 
energy dispersion curve of the single DW as 

oj{k) = \{Ji + J2)+eJicos{^k). (6) 

As was seen above, the NNN terms can not contribute to 

the low-energy dynamics of the single DW sector within 
the first order of e. Then a natural question arises: what 
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Fig. 5. Low energy dynamics of the DW excitations. The NN 
tcrm{Ji) shifts the = 1/3 DW smoothly, while the NNN(J2) 
term move the DW bound state. 
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is the role of the NNN term? In order to see it, we exam- 
ine the effect of the NNN term on the DW bound state 
of 5^ = 2/3. We define the state of the DW bound state 
at a position z by 



iTTiTTTntTi 



(7) 



Then the matrix element of TYnnn for the DW bound 
state is obtained as 

WnnnIz) = ^(^1 + J2)\^) + ^ (k - 3) + |z + 3)) , (8) 

implying that the DW bound state can move by using the 
NNN term without extra energy cost. This is a crucial 
point on the NNN term. When Ji ^ J2, the low-energy 
excitation on the plateau is basically described by the 
single DW excitation. On the other hand, when Ji <C J2, 
we can expect that the role of the DW bound state be- 
comes essential in the low-energy excitation. The switch- 
ing mechanism between the single DW excitation and the 
DW bound state is clearly associated with the crossover 
between the single chain nature and double chain nature 
of the zigzag chain. 

We further analyze the two DWs and their bound state 
problem systematically. By solving the two DW problem 
in = 2/3 sector in the bulk limit, we obtain the dis- 
persion curve of the bound state particle 



\u) = - (Ji + J2) + eJi(e"''" -(-e-^^)cos3w, (9) 



where u is the momentum of the mass center and v is 
given by e^"" — a ™^ . Since the physical solution must 
satisfy e^^ > 1 and — 1 < cos 3m < 1, the bound state dis- 
persion curve emerges in restricted ranges of u. In FIG. 6 
we illustrate the numerically obtained dispersion curve 
(9) for e = 0.1. As Ji is decreased, the bound state band 
comes down to the low energy region. In particular the 
bottom of the bound state band becomes lower than that 
of the free DW band for < 0.5. As will be discussed 
in the next section, this is important for the switching of 
the cusp and even-odd effect of the magnetization curve. 

In order to see another evidence of the = 1/3 
DW, we have calculated the dynamical structure factor 
{k,oj) for a — 0.7 and e = 0.1 at the 1/3 plateau 
with the exact diagonalization for 30 sites. Since the 
— +1 excitation is represented as a combination of 
three = 1/3 DWs, we can expect the three DW con- 
tinuum for a — 0.7. Figure 7 shows the comparison of 
the {k,uj) to the three DW continuum obtained by 



Fig. 6. The bound state dispersion of S'' = 1/3 DWs for e = 0.1 
and J2 = l.'^^ The curves are shown for < n < 2it/3. The 
sohd hnes indicate the curves for Ji =1.0, 0.8, 0.5 and 0.4. The 
dashed hnes are the lower bounds of the free two DWs. 
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Fig. 7. Dynamical structure factor {k,u)) at the 1/3 plateau 
for Ji = 1.0, J2 = 0.7 and e = 0.1. The radius of a circle is 
proportional to the intensity of the structure factor. Solid lines 
mean the upper and lower bounds of the three DW continuum. 



the 2nd order perturbation of e. Then we can see a good 
agreement, although the intensity near the lower bound 
is weak due to the finite size effect.^* 

We can also obtain the bound state dispersion curve 
of the = ~l/3 DWs below the 1/3 plateau as 



cos3m, (10) 



where e 



3v 



Q COS 6it — A 
cos 3ii 



with A 



_ -l+2a 



originating 



form the "binding energy" in the zeroth order (Ising 
limit). Here we just note that the A term yields a sensi- 
tive a dependence of the bound state dispersion around 
a = 1/2.25 

5. DW excitation and the magnetization curve 

We discuss the relevance of the DW excitations to the 
magnetization curve. The magnetization curve of the ID 
quantum spin system is generally described by the hard- 
core bosonic particle picture;'^2.33 particle having a 
magnetization fills its dispersion curve up to the "chemi- 
cal potential" or "fermi level" corresponding to the exter- 
nal magnetic field. Then the magnetization curve is in- 
terpreted as the chemical potential v.s. particle-number 
curve, where the shape of the dispersion curve is essen- 
tially important to figure out the feature of the magneti- 
zation curve. As was seen in the previous section, the ex- 
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Fig. 8. Schematic diagram (a) for the cusp singularity and (b) for 
the even-odd oscillation of the magnetization curve. The open 
circle denotes the DW excitation and the solid circle denotes 
the DW bound state. The broken lines indicate the "chemical 
potential (fermi level)" corresponding to the magnetic field 

citations around the 1/3 plateau are described well by the 
DW particles carrying = ±1/3 and their bound states 
having = ±2/3. Since the total magnetization of the 
system always takes an integer value particularly for a 
finite size system, the magnetic excitations of S'^ = ±1 
around the plateau are described by combinations of the 
fractional-S'^ DW excitations. 

We first discuss the relation between the S"^ = ±1/3 
excitation and the characteristic properties of the high- 
field branch(above the 1/3 plateau). If Ji 3> J2, the effect 
of the NNN term is not so big. Thus the low-energy ex- 
citation aromid the plateau state is basically described 
by the free S^' = 1/3 DW excitation. The magnetiza- 
tion curve near the 1 /3 plateau is reflecting the shape of 
u){k) for the = 1/3 DW; the magnetization raises from 
the 1 /3 plateau with the square-root behavior associated 
with the curvature around the bottom of cj(A;). 

As J2 is increased and the NNN effect becomes more 
significant, the dispersion CTirvc of the DW bound state 
comes down to the low-energy region. Then there are 
two possible situations: the bottom of the DW-boimd- 
state dispersion curve is slightly higher than that of the 
single DW dispersion curve(FIG.8-(a)), and the opposite 
(FIG.8-(b)). Here it should be remarked that the switch- 
ing point between the situations in FIG. 8 (a) and (b) for 
the isotropic case corresponds to a ~ 0.82, at which the 
cusp lines meet the higher field edge of the 1/3 plateau 
in FIG.l. 

For the case of FIG.8-(a) , the magnetization increases 
along the single DW dispersion curve, as long as the fermi 
level is below the bottom of the bound state dispersion 
curve. However, when the fermi level touches the bot- 
tom of the bound-state dispersion curve, the magneti- 
zation curve captures the band edge singularity to have 
the cusp(FIG.8-(a)). After the fermi level exceeds the 
bottom of the bound-state dispersion curve, the magne- 
tization can increases rapidly by using the bound-state 
dispersion curve. In the context of number of the cross- 
ing points of the fermi level and the dispersion curves, 
we can see that the one component TL liquid is realized 
between the field cusp and the 1/3 plateau, where the 
fermi level intersects only the single DW dispersion. On 
the other hand, the two component TL liquid is expected 
above the cusp singularity, since the fermi level crosses 
both of the single DW and bound state bands. 

For the case of FIG.8-(b), the magnetization curve in- 



creases by using the bound state dispersion curve. Here 
we recall that the excitation of 5^ = 1 is always repre- 
sented as a combination of three DW particles and thus 
there are 3|M - 7V1 number of = 1/3 DWs in the sys- 
tem. Then, if M =odd for N =even, 3\M-N\-1 number 
of the DWs can conform the bound state of = 2/3, 
but the remaining one DW can not find its partner. Thus 
the remaining one DW has to sit on the free DW disper- 
sion curve, as in the diagram of FIG.8-(b). Therefore the 
M =odd state has a slightly higher energy due to the 
gap between the bottom of the single DW band and the 
fermi level lying in the bound state band. On the other 
hand, for the case of M =cvcn, all of the DWs can find 
their partners and conform the bound states successfully. 
Clearly this is the origin of the even-odd behavior of the 
magnetization curve, for which the one component TL 
liquid composed of the DM bound state can be expected, 
since the fermi level crosses only the bound state band.^^ 

The analysis of the low-field branch can be done by 
almost the same line of the argument. We can easily see 
that for a <C 1/2 the shape of the bound state dispersion 
corresponds to the case (a), while a 1/2 corresponds 
to the case (b). However, in contrast to the high-field 
branch, the = —2/3 DW bound state has the zeroth 
order binding energy, namely A term, which yields the 
sensitive behaviors of the dispersion curve of the bound 
state near a = 1/2. We have calculated the magnetiza- 
tion curve of e = 0.1 around a ~ 1/2 intensively, and 
verified that the magnetization curve actually captures 
the sensitiveness of the dispersion curve. 

As seen above, the fractional-S'^ DW excitation picture 
can successfully explain the magnetization curve calcu- 
lated by the DMRG. Since the dispersion curves in §4 
are based on the first order perturbation, the higher or- 
der contribution of e may affect the shape of the disper- 
sion curves. In order to check this point, we have also 
performed DMRG calculations for various e and verified 
that the 1/3 plateau and the cusp singularities in the 
Ising limit are adiabatically connected with the isotropic 
case; the topology of the magnetic phase diagram around 
the 1/3 plateau is almost the same as the present Ising- 
like XXZ model. Thus the DW picture is also maintained 
against the quantum fluctuation of the XY term. 

Further we should make a comment on a connection 
to the other pictures of the magnetization curve. A good 
picture to describe the magnetization curve near the sat- 
uration field is the spin wave excitation from the fully 
polarized state, the shape change of which explains the 
cusp singularity successfully.-'^* Similarly, the magneti- 
zation curve near the zero magnetic field can be also 
explained by the shape change of the dispersion curves 
of the DW in the dimer-singlct states for the isotropic 
case(e = 1),2,5,37 ^j. q£ dispersion of the DW in 
the Neel type order for the sufficiently Ising anisotropic 
case."'' These excitations from the saturation field/zero- 
magnetic field limits and the fractional DW picture based 
on the 1/3 plateau are complemental to each other for a 
quantitative analysis of the magnetization curve. The re- 
lation between these pictures from the 1/3 plateau and 
the saturation/zero-magnetic fields is a remaining im- 
portant problem. In addition, we also note that how the 
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system in Ji ^ approaches the decoupled chains limit 
is also an interesting problem, since the bosonization re- 
sults may be available. 

6. summary 

We have reviewed the magnetic phase diagram of the 
S = 1/2 zigzag spin chain including the strongly frus- 
trated region and its underlying microscopic mechanism. 
Starting from the 1/3 plateau state in the Ising limit, we 
have constructed the = ±1/3 DWs and their bound 
states carrying = ±2/3. Then wc have emphasized 
the role of the single chain and double chain structures 
of the zigzag lattice; Whether the system has the single 
or double chain characters can be labeled by an integer, 
while the parameter a interpolating them is a continu- 
ous real number. Thus, as a varies gradually, the single 
and double chain natures of the system have to change 
at a critical value of a. We think that the = ±1/3 
DW excitations and their bound state strikingly demon- 
strate such switching of the the single and double chain 
natures of the excitations, which explains the character- 
istic behaviors of the cusp singularities around the 1/3 
plateau. However, the "interplay" of the integer quan- 
tum number and the continuous parameter itself should 
be independent of the microscopic details of the excita- 
tions. Wc think that this fact may be a essential view 
point behind the exotic structure of the phase diagram. 

Although the zigzag chain is a quite simple model, we 
think that it captures the intrinsic aspects common in 
a class of frustrated spin chains having the zigzag type 
structure. Recently, Z3 criticality in related models is in- 
vestigated,'^*"'*'' where the connection to the 1/3 plateau 
of the zigzag chain may be seen. The stability of the 1/3 
plateau is also an interesting problem; a perturbation 
of the trimarized interaction may induce • • • "f-singlet- 
t-singlet- • • type order rather than the tti-*^ We hope 
that our results stimulates further theoretical and exper- 
imental studies of physics induced by the frustration and 
the magnetic field. 
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